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Abstract 

We compute the contribution to the total cross section for the inclusive pro¬ 
duction of a Standard Model Higgs boson induced by two quarks with different 
flavour in the initial state. Our calculation is exact in the Higgs boson mass and the 
partonic center-of-mass energy. We describe the reduction to master integrals, the 
construction of a canonical basis, and the solution of the corresponding differential 
equations. Our analytic result contains both Harmonic Polylogarithms and iterated 
integrals with additional letters in the alphabet. 

PACS numbers: 14.80.Bn 12.38.B 



1 Introduction 


The precise determination of the properties of the recently discovered Higgs boson is 
among the main tasks of the upcoming run II of the CERN Large Hadron Collider (LHC). 
A crucial input to this enterprise is the total production cross section in gluon fusion. 

Leading order (LO) contributions to a{pp H + X) were already computed by the end 
of the 1970s in Refs, and the next-to-leading order (NLO) QCD corrections have 
been available for almost 20 years including the exact dependence on the top quark 
mass (see also Ref. for analytic results of the virtual corrections). NLO electroweak 
corrections have been computed in Ref. p(0] and mixed QCD-electroweak corrections are 
considered in Ref. 


11 


At LHC energies the NLO QCD corrections amount to 80-100% of the LO contributions 
which makes it mandatory to compute higher-order perturbative corrections. At the 
beginning of the century three groups independently evaluated the next-to-next-to-leading 
order (NNLO) corrections 12 -15 in the limit of an inhnitely heavy top quark. Finite top 

turn out to be at most 


16-22 


quark mass effects, which have been investigated in Refs, 
of the order of 1%. 

At next-to-next-to-next-to-leading order (N^LO) several groups have contributed building 
blocks to the total cross section. In Refs. [23}]^ the effective Higgs-gluon coupling has 
been computed to four-loop accuracy. In preparation of the N^LO calculations the 0(e) 
contributions to the NNLO master integrals have been computed in Refs. 22 ,^ where 
d = 4 — 2e is the number of space-time dimensions in dimensional regularization. Results 
for the LO, NLO and NNLO partonic cross sections expanded up to order e^, and e^. 


respectively, have been published in Refs. 27,28 . All contributions from convolutions 


of partonic cross sections with splitting functions, which are needed for the complete 
N^LO calculation, are provided in Refs. 


expression has been constructed in Ref. 
for Us 


27-29 . The full scale-dependence of the N^LO 


Three-loop ultraviolet counterterms needed 
30l 31 and the operator in the effective Lagrangian were computed long ago. 


Within the effective theory, three-loop virtual corrections to the Higgs-gluon form factor 
have been obtained by two independent calculations 


33,34 


single-soft current to two-loop order has been computed in Refs. 36,37 


(see also Ref. [^). The 
which is an 

important ingredient to the two-loop corrections with one additional real radiation. The 
latter have been computed in Refs. 38,39 . The single-real radiation contribution which 


originates from the square of one-loop amplitudes has been computed exactly in terms 
of the Higgs boson mass and the partonic center-of-mass energy in Refs. 40,41 . The 


soft limit of the phase space integrals for Higgs boson production in association with two 
soft partons were computed in Refs. 42,43 , in the latter reference even to all orders 


in e. The triple-real contribution to the gluon-induced partonic cross section has been 
considered in Ref. 44 in the soft limit. In particular, a method has been developed which 
allows the expansion around the soft limit. A similar analysis for the double-real-virtual 


contributions has been published in Ref. 45 
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The two leading terms in the threshold expansion for the complete N^LO total Higgs pro¬ 
duction cross section through gluon fusion have been presented in Refs. 42,46,47 . How¬ 


ever, for physical applications more terms in the threshold expansion are necessary 
In fact, in Ref. 


46 


48 more than 30 expansion terms have been computed which is sufficient 


for all phenomenological applications. It is important to cross-check the result of Ref. 48 


In this paper we present the first step in this direction. In particular, results are obtained 
which are exact in the Higgs boson mass and the partonic center-of-mass energy. 

Further activities concern the development of systematic approaches to compute the mas¬ 
ter integrals for a{pp H + X), see, e.g.. Refs. |3^[40|[4l|[44|[4^ . 

Several groups have constructed approximate N^LO results for the total cross section 
taking into account information from the soft-gluon approximation and the high-energy 


limit 50 56 


In the following, we briefly outline the framework which we use for our calculation. In 
the limit of an infinitely heavy top quark the effective interaction of the Higgs boson with 
gluons is described by the Lagrange density 


C 


H 


y,efr 


,( 5 ) 




( 1 ) 


where is the usual QCD Lagrange density with five massless quarks, H denotes the 

Higgs field, V its vacuum expectation value and Ci is the matching coefficient between 
the full and the effective theory. is the gluonic field strength tensor constructed 

from fields and couplings already present in The superscript “0” denotes the bare 

quantities. Note that the counterterms of H^/v^ are of higher order in the electroweak 
coupling constants. 

The top quark mass enters the cross section via the matching coefficient Ci whereas the 
quantities in the effective theory depend on 


mt 


X = 


( 2 ) 


where rrih is the Higgs boson mass and \/§ the partonic center-of-mass energy. For later 
convenience we also introduce the variable 


y = 1 — X . 


( 3 ) 


At the partonic level several sub-processes initiated by quarks and gluons in the initial 
state have to be considered. The numerically most important but also technically most 
complicated contribution is the one with two gluons in the initial state. In the present 
paper we consider the subprocess qq' —)■ H + X at NNLO and N^LO. Its phenomenological 
impact is very small, but we use this process to demonstrate our method which leads to 
exact results in x and avoids the high-order soft expansion. 
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Figure 1: Sample Feynman diagrams for qq' —>■ qq'. The imaginary parts due to Higgs 
boson cuts provide the cross section for the process qq' ^ H + X at NNLO and N^LO. 
Solid, curly and dashed lines represent quarks, gluons and Higgs bosons, respectively and 
blobs denote the effective Higgs-gluon couplings. 


For the calculation of the total cross section it is convenient to consider the imaginary part 
of the forward scattering amplitude qq' —)■ qq'. Sample Feynman diagrams contributing at 
NNLO and N^LO are shown in Fig. [Tj To obtain the cross section all cuts involving the 
Higgs boson have to be computed which means that both three- and four-particle cuts 
have to be considered at N^LO. There are no two-particle cuts. 

The remainder of the paper is organized as follows: In the next Section we discuss the re¬ 
duction of the full set of integrals to master integrals and the construction of the canonical 
basis. For the latter integrals a system of differential equations is derived. The following 
two sections are dedicated to the evaluation of the initial conditions involving cuts of three 
(Section and four (Section particles. In Section we introduce recursively dehned 
iterated integrals which are needed for the analytic representation of the hnal result. The 
partonic cross section is discussed in Section where analytic results are given. Finally 
we conclude in Section 0 


2 Reduction and canonical master integrals 


We generate all two- and three-loop forward-scattering amplitudes for the process 
q{Pi)q'{P 2 ) —^ q{Pi)q'{P 2 ) involving a virtual Higgs boson with the help of qgraf 
and process the output hie to select the contributions which contain cuts through the 
Higgs boson line. This leads to 1 two-loop and 224 three-loop Feynman diagrams. At 
three-loop order the corresponding amplitudes can be mapped to 17 integral families 
which are shown in Fig. For each of them reduction tables are generated using a com¬ 
bination of the publicly available program FIRE 58 and in-house programs, rows and 
TopoID 59 , which implement the Laporta algorithm 1^. The use of rows and TopoID 


guarantees that all available symmetries are exploited which is important to minimize the 
number of master integrals. After completing the reduction for each family we obtain 332 
master integrals. In our next step we minimize the number of integrals by simultaneously 
considering all families which leaves us with 111 master integrals, 108 of which are needed 
for the cross section. In the following we refer to this set of master integrals as “Laporta 
master integrals”. 


4 
















10 8 



7 4 

BTl 

8 4 


2 

1 

6 


9 

10 


BT5 

6 7 


3 

1 

8 

2 

5 



4 


10 

9 



BT9 

9 



BT13 


7 4 



9 8 

Q 


10 6 
BT6 


10 5 



7 



8 5 




10 



BTll 


9 



9 7 


5 8 

BT4 

9 10 


8 5 
BT8 

9 


8 6 
BT12 


8 







9 



Figure 2: Graphical representations of the 17 three-loop integral families. Plain and 
double lines indicate massless propagators and the Higgs boson lines, respectively, and 
the wavy lines indicate the possible cuts. 


Note that we have performed the calculation for general gauge parameter ^ which drops 
out after relating master integrals from the different families. This constitutes a strong 
check on the correctness of our result. 

For the evaluation of the master integrals we follow the ideas of Ref. and construct a 


5 










































































canonical basis which allows for a simple and straightforward solntion of the corresponding 
differential eqnations (see Refs. for reviews on the nse of differential eqnations for 

the compntation of Feynman integrals). Whereas most of onr calcnlation is antomated 
to a high degree the constrnction of the canonical basis reqnires mannal manipnlations of 


each individnal integral. We have applied several tricks described in the literatnre 64 68 
and also follow the algorithm developed in Ref. 
canonical master integrals in conpled snbsystems. In Ref. 69 


which allows the constrnction of 
an algorithm has been 
snggested which antomates the constrnction of the canonical basis. However, a pnblic 
implementation is not yet available. 


In a canonical basis the differential eqnations have the form 


dxf(,x,e) = eA{x)f{x,e), 


( 4 ) 


where f{x,e) is a vector containing all canonical master integrals. In onr case the matrix 
A[x) can be written as 


A{x) 


abed e 

- -h -1-h - \ - , , 

1- X 1+ X X 1+4:X x^/l + 4:X 


( 5 ) 


where a,..., e are constant matrices. The hrst three terms on the right-hand side of 
Eq. ([^ lead to the well-known Harmonic Polylogarithms (HPLs) fTOl (see Refs. fTTlI^ 


for a convenient Mathematica implementation) in the solntion of the master integrals. 
The fonrth and hfth terms in Eq. ([^ are only needed for the integral family BT3 as we 
will describe in detail in Section HI 


Besides the simple solntion of the differential eqnations the canonical basis also has the 
advantage that for the initial conditions only the leading terms of order are needed in 
the soft limit. As a conseqnence, no explicit calcnlation is needed in case the hrst non-zero 
contribntion of a canonical master integral is of 0{y) or higher. In onr calcnlation the 
bonndary conditions are compnted for Laporta master integrals. Afterwards the resnlts 
are transformed to the canonical basis. 


3 Three-particle cuts 


The three-particle-cnt contribntions contain a one-loop snb-diagram. As onr hrst step 
we represent the loop in terms of Mellin-Barnes integrals and perform the momentnm 
integration. Afterwards in the soft limit all integrals are represented as phase space inte¬ 
grals of soft partons, which can be converted to integrals over energies and angles. These 
integrals are also calcnlable nsing Mellin-Barnes integrals. Hence, we obtain mnltifold 
Mellin-Barnes integral representations for each master integral in the soft limit. They are 
evalnated extending the method developed in Ref. 


73 


for the calcnlation of the three- 
where fonr-particle 


44 


loop static potential. The notation is mainly adopted from Ref. 
ents have been considered. In this reference also a techniqne has been developed which 
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transforms soft phase-space integrals to Mellin-Barnes integrals, which has been applied 


in Ref. to three-particle phase space contribntions. In contrast to Ref. 45 we do not 
apply the method of regions to compnte the integrals. 

Before describing the procednre in more detail we have to introdnce some notation. We 
denote the external momenta by Pi,P 2 and the light momenta involved in the cnt by 
P 3 ,P 4 and p 5 , where p^ will occur for the four-particle phase space integrals in Section]^ 
Loop momenta are denoted by Vi. According to Ref. 44 the scaling of the phase space 

y^/s for i = 3,4 and 5 


momenta in the soft limit is given by pi ~ ^/s for i = 1, 2 and pi 
in the center-of-mass frame of the incoming quarks. We eliminate the momentum of the 
Higgs boson in favour of the momenta of the massless partons and dehne rescaled scalar 
products 


Sij — 


sy ' 
{pi+pjf 

sy^ 


z = 1, 2 and j > 2 , 
i > 2 and j > 2. 


( 6 ) 


Furthermore, we use the energies and angles parametrization 

\T 


ZL 

\fs 

Zz 

Pi 

Vsy 


2^^! “ 1) 

2/^2 = -(l,0d_2,-l)'''. 


\ErPi 


for z > 2 , 


( 7 ) 


where Ei parametrize the partons’ energies and their d-dimensional velocities. 0^-2 is 
an abbreviation for a sequence of d — 2 zeros. For later convenience we also introduce 
f^ij 

In the following we exemplify the individual steps of the algorithm on the integral 
Bg = BT9(1,1,1, 0,1, 0,1, 0,1,1,0,0) 


d^>3 


N 


{27r)‘^v‘^{Pi-P3y{Pi-P3-Pi + vy{Pi + P 2 -P 3 -P 4 + ■ 

is a normalization factor given by 


N = — 


(dvr; 


i2-e \ 3 


( 8 ) 


( 9 ) 


27r \F(1 -|- e) 

where the factors F(1 -|- e) and {dnY are introduced for convenience and d^g is the soft 
three-particle phase space measure which can be written as 


d<f>3 = (27r) 


-5+4eo-6+4e„l-2e„,3-4e 




(i-t^.)riF-/ 


dEi I dllf-'. (10) 


ilf ^ is the d-dimensional solid angle. 

The algorithm for the computation of the three-particle cut contribution is as follows: 


7 









1. Introduce a regularization parameter S for the numerators. This is necessary to 
avoid terms r(0) which otherwise could appear in step below. We introduce 6 to 

2 r 2 ”1 

the exponent of the scalar products, namely, {pi + ■ ■ ■) hm^^o [{Pi + ''') J 

2. Perform sub loop integration and introduce Mellin-Barnes integrals. 

We (i) introduce Feynman parameters to combine propagators involving loop mo¬ 
menta, (ii) perform loop integration and (iii) introduce Mellin-Barnes variables to 


obtain a factorization of the Feynman variables 44 using the formula 

r-+ioo 


1 


(X + YY F(A) 27rz 


d^F(A + ;2)F(-;2) 






( 11 ) 


In our example we obtain a one-fold Mellin-Barnes integral over zi which has the 
following form 

o7-4e.^3-2e 


_ 7r=^-^T(-^i)F(^i + 1)F(-6 )F(-^i - 6)F(^i -F 6 + 1) 

(Pl-P3)^(Pl +P2-P3-P4)"^^i(Pi-P3-P4)^^^+^"+^r(l - 2e)F3(e -F 1) 


( 12 ) 


3. Express the propagators in terms of velocities and energies, and take the soft limit, 
i.e., y -hO. 

Using Eq. we can replace the propagators in our examples Bg as 


(Pi-Ps^ 




syE3/313 ’ 


(P,_P3_)j2+2«+2. ^ (-ls!/BsA3-ii>!/B4ft4 + iVBsB4/334 

(p, +p,_)3_p,)- 2.. ^ (yA2 - i3Si53A3 - ii>!/B4A4 - ^i>!/B3A3 

1 1 A""' 

--syEilSgA + -sy'^E3E4/334j ■ 

To leading order in y this becomes 

1 2 
-)■ 


—zi—e —1 


{pi-pYj"^ 


syE3l3i3 ’ 


1 


-)■ 


—syE3j3i3 - ]^syEif3iA 




zi 


( 13 ) 
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4. Introduce Mellin-Barnes variables to factor the and E dependence. 

In our example a further Mellin-Barnes parameter Z 2 has to be introduced to de¬ 


compose the sum on the r.h.s of Eq. (14). Afterwards, the energy integrations are 
trivial and we obtain 


Bo = 


r(l-2e)r3(e + l) J 2TriJ 2nV ’ ^ ^ ^ 


X 


r(— Zi — 6 )r(x 2 — 2e + 2)r(— Zi — Z 2 — 3e)r(xi + X 2 + e + 1) 

r( —Zi — 5e + 2) 




(15) 


5. Convert angular integrations to Mellin-Barnes integrations. This is achieved by 
using repeatedly 


da" 


d-i 




Jni 

2^“Em=l «m —2e.^l—e 


nLi r («.) r (2 - e::.=i - 2e) 


n n n 

——T ( 


X 


dx 




2 fci 

jkjl 


_fe=l l=k 

in order to perform the fl integrations. 
In the case of our example this leads to 


nr 

.k=l 


m=l 


(Xp, 


k=l l=k 


+ ^ Zkl 

l=k 


1=1 


,( 16 ) 


Bo = s 


-3e-l 


H ^ /* H ^ 

^ + i)r(- 2 ; 2 ) 


X 


r(—e)r (—Zi — e)r(x 2 — e -l- l)r(— Zi — Z 2 — 2 e — l)r(xi -|- X 2 + e -|- 1 ) 


r(l - 2e)r3(e + l)r(-xi - 5e + 2) 


(17) 


6 . Simplihcation with Barnes’ Lemma. 

We use the routine DoAllBarnes [] of the package barnesroutines .m [^. Before 
applying it, we convert the cosine to Gamma functions using either 

cos(a) = —7 -- 7 - 

r(i-f)r(j) r(i-j)rq) 

or 

71 

cos(a) = -T-j-r- 71 -^ (19) 

r(i-f)r(i + f) 

depending on whether half-integer arguments are present in the hnal expression or 
not. The latter should be avoided to arrive at simpler expressions. 
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7. Take the limit 5 —?• 0 (if needed) and expand in y and e. 

Using MBcontinuef] from the package MB.m |[^, we can obtain Mellin-Barnes rep¬ 
resentations for the limits 5 —)■ 0 and e —)■ 0. To achieve this goal, we have slightly 
modihed the code to prevent that log(j/) terms appear. 


After that we expand the representation in 5 and e using MBexpand [], and in y 


using MBasymptotics [] 77 


8. Further simplihcation of Mellin-Barnes integrals. 
We apply the following procedures iteratively: 


• MBapplyBarnes[] 

• DoAllBarnes [] 

• Simplihcation of the integration contours such that all integrals with the same 
number of Mellin-Barnes parameters have the same integration contours. 


9. Conversion to nested sums and their evaluation. 


To achieve this, we hrst use the residue theorem to convert the integrals to sums. 
In case this step generates divergent inhnite sums, we introduce a regulator 
in the integrand, where the q’s are properly chosen numbers, cr is a regularization 
parameter, and the ZiS are Mellin-Barnes parameters in the expression. For the 
evaluation of the sums, we use the summation program described in Ref. 


73 


The hnal result for the integral Rg reads 






1„,2-56 I l„,2-4e 


+ 


15 2-5e 
4 y 


+ 


2-4e 


+ 


y 


2-,, (IK. _ 115) + ,14 _ 


V 


2-5e 


A65C2 


y 


V 4 

2-5. A925C2 


+ 18C3 


1875 \ 


16 


- 


+ y‘-*‘ (- 36(2 - 11(3 + 60 ) 


V 8 


+ 135(3- 


8 


32 


+y 


2-4e 


-I68C2 - 66C3 + + 248 


+ e 


y 


2-5. 


-I98C3C2 + 


20625(2 ^ 1575(3 465(4 


+ 294(5 


196875 \ 


16 ' 2 16 ' " 64 y 

+ y^-^^ ( 132 ( 3(2 - 720(2 - 308(3 + 315(4 - 105(5 + 1008) 

0{e^) + 0{y ^), (20) 


where terms up to 0{e^) have been computed. For brevity only terms up to order are 
shown. 
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We have used the described algorithm for all needed three-particle initial conditions with 
one exception: the result of the integral BT9(1,0,1, 0,1,1,1, 0,1,1, 0, 0) where the lines 


{1, 7, 9} are cut is taken over from Eq. (5.32) of Ref. 45 


As a cross check we have computed more integrals in the soft limit than actually necessary 
to £x the boundary conditions. Afterwards we have checked that the solution of the 
differential equation reproduces these additional terms. 


Note that the algorithm described in this section can also be applied to the four-particle- 
cut contribution after applying obvious modihcations. In this way we have cross checked 
most of our results, which have been obtained using the method which we describe in the 
next section. 


4 Four-particle cuts 


To compute the initial condition of the four-particle-cut contributions we closely follow 
the procedure described in Ref. (^. For completeness we briefly repeat the individual 
steps in this section. The soft expansion of the four-particle cut integrals exhibit only one 
region, which is dehned by the scaling in y of the scalar products Sij dehned in Eq. (|^. 
Reversed unitarity allows for an expansion in the limit y —>■ 0 of the Higgs boson 
propagator which in our parametrization is given by 


(Pl +P2-P3-P4- PbY - X 

OO ^ 

= ~ [“('^34 + 535 -f 545)]^ 

fc =0 ^ 


(1 + S13 -|- S23 + Si 4 -|- S24 + -Sis -|- S25) 


k+1 


( 21 ) 


where the subscript “c” reminds that the propagator has to be cut. In the soft limit 
only the term /c = 0 is needed. The massless propagators of the quarks and gluons are 
expanded as a Taylor series in the limit r/ —)■ 0 as well. This yields shifts in indices of 
the propagators, which are removed by subsequently applying the Laporta algorithm 
as implemented in FIRE 58 in the soft kinematics. We obtain eleven master integrals. 
Ten are given in Ref. 44 where analytical results are derived. The eleventh integral 


corresponds to the soft limit of BT1(1,1,1, 0, 0, 0,1,1,0,1, 0, 0) (cf. Fig. which can be 
cast in the form 


’ll 




dd)^ 


(Sl3 -f Sl4)(si4 + S 15 ) 


( 22 ) 


where <h 4 is dehned in analogy to in Eq. (10). In Ref. this integral probably only 
contributes to higher orders in y which is why it has not been discussed in that paper. 


Following Ref. 44 we apply Eq. (11) to convert the sums in the denominator of Eq. (22) 


into products at the cost of introducing Mellin-Barnes integrals. 


11 














Introducing energies and angles in analogy to Eqs. (§ and ([^ one can integrate the 
energies in terms of F functions, such that the only non-trivial integrations are given 


by three integrations over solid-angles, each of the form of Eq. (16), which are turned 


into Mellin-Barnes integrals. Following this procedure, we arrive at a one-dimensional 
Mellin-Barnes integral 


^+200 


Fn(e) = 


dz r(6 — 6e)r(—2e — z)T{—e — 2 :)r(—z)r(l -(- 2 :)r(l — e + z) 


2'Ki 


r(4-6e)r2 (l-e)r(l-2e-.2) 


.(23) 


which we expand in e and solve by applying the algorithm of Ref. As a hnal result 
we obtain 


Fii(e) = 20 C 2 + e(-54C2 +140(3) +6^(36(2 - 378(3 + 600(4) 

+ 6^(252(3 + 160(2(3 - 1620(4 + 1860(5) 

+ €^(-432(2(3 + 560(1 + 1080(4 - 5022(5 + 6420(6) 

+ 6^(288(2(3 - 1512(1 + 4800(3(4 + 3348(5 + 960(2(5 - 17334(6 + 15240(7) 
+ 0 { e ^), (24) 


which we have checked numerically using the package MB.m 76 
the integral^ + 2 ( 6 ),..., + 10 ( 6 ) of Ref. 44 

cients of Eq. ( [24l ) are integers, an observation also made in Ref. 

+2(6),..., +10(6). 


We have also rederived 
It is interesting to note, that all coeffi- 
|44) for the integrals 


For many master integrals, we computed more terms in the soft expansion than required 
to hx the integration constants. These terms could be compared to the expansion of the 
exact result and thus strong consistency checks are obtained. 


Note that an alternative method to compute four-particle phase-space integrals in the soft 
limit has been developed in Ref. [78]. 


5 Iterated integrals beyond HPLs 


The solution of 16 out of our 17 families can be expressed in terms of HPLs 70 , however. 


for BT3 this is not possible. In fact, the differential equation of the canonical basis implies 
an alphabet for the iterated integrals which involves square roots. The letters are 


(25) 


x’l — a:’l + a:’l + 4a:’ xy/l + dx 


The master integrals in which the last two letters show up can be classihed as having the 
common subtopology drawn in Fig. The contributing integrals with this property are 

BT3(1, 0,0,1,1,1,1,1,1, 0,0,0), 


Whe integral J^i(e) is simply the volume of four-particle phase space itself. 
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BT3(1, 0,-1,1,1,1,1,1,1, 0,0,0), 

BT3(1, 0,0,1,1,1,1,1,1,1,0,0), 

BT3(1,1,1,1,1,1,1,1,1, 0,0,0), 

BT3(1,1,1,1,1,1,1,1,1,0,-1,0), 

BT3(1,1,1,1,1,1,1,1,1,-1,0,0). (26) 


In general the occurrence of square roots yjx — Xq can be anticipated by observing half¬ 
integers in the diagonalized form of the matrix residue at Xq, as shall be briefly explained in 
the following using the above example. Let us denote the system of differential equations 
for the integrals in Eq. (26) by dxf = Af. We expand the matrix elements of A in a 
Laurent series around xq = —1/4 and take the coefficient of {x — xq)”^, which is called 
the matrix residue. After diagonalization we obtain 


diag(0,0,0,0,0,l/2-4e). 


(27) 


In a next step we expand the element /g corresponding to the last entry of the diagonal 
matrix in a power series in {x — Xq) and obtain with the help of the differential equation 
dxfe = (1/2 — 4e)/6/(x — xq). Note that the occurrence of the half-integer prefactor on 
the right-hand side (for e —?• 0) implies the occurrence of the square root \/l + 4x, which 
in the full solution may show up in coefficients and in the alphabet of iterated integrals. 


The present calculation contrasts earlier ones encountering square root letters 79-81 


where the occurrence of the square root is connected with the presence of massive two- 
particle or four-particle cuts in the integrals (cf. the connection of square root letters in 
iterated integrals with (inverse) binomial sums in 
man diagrams involving (inverse) binomial sums in 


82 


as well as calculations of Feyn- 
Topology BT3, however. 


83-86 


represents four-particle phase-space integrals with only one massive line in the hnal state. 



Figure 3: Common subtopology of all the graphs in BT3 which generate square root 
letters. 

The canonical differential equation can be solved, as usual, order-by-order in e. Afterwards 
the constants of integration have to be determined. This is done by expanding the generic 
solution in a generalized Taylor series expansion around x = 1 and matching with a 
calculation in the soft limit. For the expansion one needs to extract the logarithmic part 
due to log(l — x). This can be done using the shuffle algebra, and making sure that 
1/(1 — x) never occurs in the rightmost index of the iterated integrals. As a result the 


13 




















iterated integrals either diverge like log(l — x) or are regular in the limit x 1. For the 
matching procedure one now only needs the (1 — x)°-order, i.e. the regular part evaluated 
at a: = 1, while logarithmic orders provide a cross check for the generic solution with the 
calculation of the boundary conditions. 


In this way the canonical master integrals and hence the Laporta masters are expressed 
in terms of iterated integrals over the alphabet (25). For numerical evaluations, it is 
advantageous to modify the above alphabet to be 


111 11 
x’ 1—x’ l + x’ l + dx’a; 

so only one letter is singular as a; —)■ 0. 

The contributions to the single Higgs boson production amplitude do not span the full 
space of functions generated by the above alphabet. In fact the relevant iterated integrals 
involving the square root letter can be constructed from 


\/l + 4x 


- 1 


(28) 


fo = 


1 

X 


f-1 


1 

1 + X ’ 


fs4 


1 

X 


( ^ 
Va/ 1 + 4x 



(29) 


For the treatment of algebraic relations and the series expansions of the iterated integrals 
with square root letters, the package HarmonicSums was used 


87-89 


For the numerical 

implementation, the convergent series expansions around x = 0 and x = 1 are helpful, 
which are available once the letter 1/x is shuffled away from the rightmost position in the 


indices of the iterated integrals. Unfortunately in contrast to the case of HPLs 70,90 


the series expansion around x = 0 has a radius of convergence of 1/4, thus more terms in 
the expansions are needed. 


The iterated integrals involving square root letters were implemented numerically in 
Mathematica, using series expansions for functions of weight <3 and up to twofold nu¬ 
merical integrals. In this way we are able to yield 10 good digits at the timescale of a 
second and below for the most complicated functions at weight 5. 


6 Results 


The total partonic cross section can be written as 

a = Cfa , 


(30) 


where Cf and a are separately finite after renormalization and the convolution of the 


lower-order cross sections with the splitting functions 27|-|29|. 
Our final result can be cast in the form 


a{qq'^H + X) = A 


a. 


TT 


42 ) 


— + — 


qq 


a. 


TT 


a 


( 3 ) 

qq' 


(31) 
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where A = Gi7’7r/(32-\/2) and is given in Eq. (54) of Ref. 14 and reads after 
identifying renormalization and factorization scale with the Higgsooson mass (i.e. /ir = 
fif = ruh) 


~( 3 ) 


+ 2fH^{x) + ^ (llx^ - 40x + 188) H^{x) - + 2fC2H^{x) 


243 


27 


1 16 

- — (772x2 + 1156x - 2213) H^{x) - — (x - l)(x + 3)H,{x)H^{x) 

24:0 9 

1 40 

- — (3139x2 + 5218x - 7620) Ho{x) - —(x - l)(x + 3)Hf{x)Ho{x) 

+ 7 ^ (83x2 264 x + 148) i7o, 1 , 1 ( 2 :) + ^ (107x2 + 344x + 164) Ho^i{x)Ho{x) 

27 ol 

32 70 

- ^ ( 2 :^ + 2 : + 4 ) C 2 Ho{x) - - l)(13x + 33)Hi{x)Ho{x) 

- ^(x + 2)2i7o,o,i(2:)i7o(2:) + ^(x + 2)2i7o,i,i(2:)i7o(2:) - ^(x + 2 ) 2^32 

—^( 2 : + 2Y(3Hq{x) —^(x — l)(x + 3)i7f (x) — 77 ^( 2 : — l)(2549x + 15343) 


27 

1 


27 


324 


11 


-(x - l)(319x + 771)i72(a;) - -(x + 2yH^^,{x) - -(x + 2)2^770,1(2:) 


27' 


54 

+ ^ (215x2 - 904x - 1402) (2 + 7 ^ (469x2 + 1840x - 218) 77o,i(x) 

- ^^( 2 : - l)(169x + 519)77i(x) - ^ (119x2 + 464x + 260) 77o, 0 , 1 ( 2 :) 

56 22 16 

+ ^( 2 : + 2 ) 2770 , 0 , 0 , 1 ( 2 :) + + 2)^77o,i,o,i(x) + ^( 2 : — l)(x + 3)C277i(x) 

+ ^( 2 : + 2)277o,i(x)77g(x) + -(x + 2)277o,i,i,i(x) + — (l3x2 + 184x + 4) (3 
9 9 oi 

1 112 

H-(l064x^ — 2853x2 + 107433x — 41149) 77 q(x) -—(x + 2)277 o,i(x)77q (x) 

729 27 

- ^ (13x2 + 135x + 164) 77o"(x) + ^ (29x2 - 28x + 28) C277o^(2:) 

+ Jr(x + l)(97x - 294)77_i(x)77o(x) + ^( 2 : - l)(97x + 285)77i(x)77o(x) 

ol ol 

- ^ (175x2 - 308x + 216) H^{x) “ “ 2)^77o,-i(x)773(x) 

- ^( 2 : - 3)(x + l)77!i(x)77o2(a;) + ^(x - l)(x + 3)HI{x)HI{x) 


27 

1 


27 


, (46480x=^ - 286656x2 + 3753336x - 1756017) 772 (x) 

OOO.Z 

35 2 

+ 7 r—(118x + 85)\/4x + 177 o2(x) - — (105x2 - 1548x - 1900) C277o2(x) 
216 


81 


15 







- ^ {9x^ -2x-8) Hsi{x)HQ{x) + —(x - 3)(x + l)H_i^s 4 (x)B'o(x) 

54 y 

35 1 

- {x — 6 )xHqs 4 { 1 )Hq{x) -(l083a;^ — 2556a; — 3046) Hq_i{x)Hq{x) 

9 ’ 81 

+ — (a; + 1) (80a;2 + 12889a; - 46117) H_i{x)H^{x) 

35 1 

——(a; — 2 )‘^Ho_i^s 4 {x)Hq{x) + -^{x — 1) (400a;^ — 1019a; + 35945) Hi{x)Hq{x) 
9 54 

35 1 

——(a; — 6 )xHq^s 4 {x)Hq{x) + — (l220a;^ — 4844a; — 7485) Hq^i{x)Hq{x) 

9 ’ 27 

+ '^{x - 2fHo_i_i{x)Hl{x) + ^{x - 3)(x + 1)770,^(1)77-1 (a;)i7o(a;) 

2 < 9 

- ^ (317a;^ - 1612a; + 212) Hofl^i{x)H^{x) “ + 2)^77o,i,i(a;)i72(a;) 

8 592 

+ - (l3a;^ + 31a; + 30) C 3 Hq{x) + —{x - 3)(a; + l)H^^{x)Ho{x) 

9 oi 

3212 4 

+ - l)(^ + 3)H^{x)Ho{x) - — (361a;^ - 1044a; + 36) CIHq{x) 

+ ^ log'^(2)(a; - 2)^77o(a;) - ^(a; + 1) (l04a;^ + 1633a; - 7915) i7!i(a;)i7o(a;) 
224 /1\ 14 

+ ( - j (a; - 2fHo{x) + ^(^ + 1) (280a;^ - 1213a; - 5327) i7_i(a;)i7o(a;) 


+ -^{x — 1) (592a;^ — 847a; + 53333) H\{x)Hq{x) —~ 2)^77g _;^(a;)i7o(a;) 

- (421040a;3 - 394707a;2 - 12461502a; + 7407221) 77o(a;) 

1 56 

- - (79a;2 - 1028a; - 212) Hl^{x)H^{x) - - \og\2){x - 2)%Ho{x) 


- — (3008a;^ - 10503a;^ + 202176a; - 46836) C 2 Ho{x) 

35 _ 35 

+ —V4a; + l(118a; + 85)i7^4(a;)7ro(a;) - — (9a;^ - 2x - 8) Ho,s 4 {^)Ho{x) 

Q TCi 

- -{x - 2)2(2770,.4(l)77o(a;) + y (a: - 2)277_i,,4(l)77o,.4(l)77o(a;) 

140 140 

+ - 3)(a^ + 1)77-1,M,s4(a;)77o(a;)-^(a; - 3)(x + l)77o,-i,M(a;)77o(a;) 

70 70 

+ — log(2)(a; - 2)277o,o,s4(l)77o(a;) “ y “ 6)a;77o,o,s4(l)77o(a;) 

70 70 

+ — (a; - 6)xHofl,s4(x)Ho(x) - — log(2)(a; - 2)^770,s4,s4(l)^o(a^) 

9 9 

70 70 

+ — (a; - 6)a;77o,M,s4(l)77o(a;) - -^{x - 6)a;77o,s4,s4(a;)77o(a;) 

9 9 

70 70 

- -(a; - 2)277-1,o,o,s4(l)77o(a;) - y (a; - 2)277-i,o,M,s4(l)77o(a;) 


16 



70 70 

- -(x - - -{x - 2)^770,-i,s4,M(l)i^o(a:) 

70 70 

—^( 2 ^ ~ Hq _i^s4^^s4:{x)Hq{x) + —(a^ — 2)^77o^o,-i,s4(l)-f^o(a^) 

y y 

H— ~ HQfl _i{x)Hl{x) — —{x — 1) (908a;^ + 7916a; — 135025) Hi{x)Hq{x) 
9 ’ ’ 81 

35 4 

- ^ (9a;^ -2x-8) Hs 4 {xYHo{x) - + l)(373a; - 1101)C2-f^-i(a;)i7o(a;) 


- — (306a;^ - 474a; - 497) Hq_i_i{x)Hq{x) - —(x - l)(307a; + 963)C2-f^i(a;)-f7o(a^) 

O -L ^ i 

1 35 

+ — (I6a;=* - 5013a;2 + 4140a; + 28873) 77o,_i(a;)77o(a;) -—(x - 2Y(2 HomMY 

z4o y 

740 8 

+ ~ '^YC 2 Ho-i{x)Hq{x) + -(a; + l)(27a; - 79)77_i(a;)77o _i(a;)77o(a;) 

70 8 

- -(a; - 2YHoAYHoM^)HYx) “ “ l)(25a; + 117)77i(a;)77o,_i(a;)77o(a;) 

+ ^ (2960a;3 + 18009a;2 - 2034a; + 69166) 77o,i(a;)77o(a;) 

448 8 

- —(a; - 2)^77o _i(a;)77o,i(a;)77o(a;) + g + l)(37a; - 113)77_i(a;)77o,i(a;)77o(a;) 

628 n 2 

+ + ‘4^C2Hq^i{x)Ho{x) + —{x - l)(475a; + 1383)77i(a;)77o,i(a;)77o(a;) 

140 16 

+ ~ 2)^77o,o- i,s4(a;)i7o(a;) - — (43a;^ - 160a; - 111) 77o _i,i(a;)77o(a;) 

496 4 

- “ 2)^77o,o- 1 -i(a;)77o(a;) + — (81a;^ - 620a; - 112) 77o,o-i(a;)77o(a;) 

896 2 

+ ~ 2)^77o,o- i,i(a;)77o(a;) - — (3996a;^ - 17088a; - 31681) 77o,o,i(a;)77o(a;) 


81 

16 


- 52(a; - 2)^77o,o,o-i(a;)i7o(a;) - — (43a;^ - 160a; - 111) 77o,i _i(a;)77o(a;) 

2 592 

+ ^ (923a;^ - 8188a; - 6336) 77o,i,i(a;)77o(a;) - - 2)^77o- 1-1 _i(a;)77o(a;) 

4 896 

+ - (275a;^ - 820a; + 484) 77o,o,o,i(2^)^o(a^) + “ 2)^77o,o,i-i(a;)77o(a;) 

64 4 

+ — (9a;^ - 20a; + 36) 77o,i,o-i(a;)i7o(a;) + — (37a;^ - 838a; + 842) C,^Hq{x) 


+ - (a;^ - 1340a; - 524) 77o,i,o,i(a;)77o(a;) 


3212 

~w' 


X + 2Y Hqi^ii{x)Hq{x) 


9 

196 500 

+ — log(2)(a; - 2 YCMx) - - l)(a; + 3)C2i^?(a;) 

+ ^(a; - 1) (256a;2 - 237a; + 23559) 773(a;) + ^ (2269a;2 - 1338a; + 9038) C; 

35 35 872 

+ (a; _ 6)a;77o,.4(l)" + 77(^ - 6)a;i7o,.4(a;)' + —(a; - 3)(a; + l)C 2 i^!i(a;) 

y y z ( 
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1208 1 

+ - l)(a; + ?>)Hl{x) - - 1) {2U48x‘^ + 200029a; - 3186149) H'^{x) 

ol 291o 

+ — (63a;2 - 6a; + 214) H^_^{x) - (494x^ - 5212a; - 4845) H^i(x) 

81 ’ 27 ’ 

1 112 

+ - 1) (57136a;2 - 1139639a; + 29021665) + — log2(2)(a; - 3)(a; + 1)C2 

26244 27 

1 56 

- (8032a;^ - 239799a;2 + 1426872a; - 432084) C 2 - ^ log^2)(a; - 3)(a; + 1) 

1458 81 

- (a; - 3)(x + 1) + ^V4a; + l(118a; + 85)Ifo,s4W 


+ ^(x-3)(x + 1)C277 o,s 4(1) - p (9a;2 - 2a; - 8) II,,(x)IIo,s4W 

140 140 

-^(a; - 3)(a; + l)i7_i,M(l)i7o,M(l) + - 3)(a; + l)i7_i,M(a;)i7o,s4(l) 

35 _ 35 

- — V4a; + l(118a; + 85)IIo,s4(^) + 7 ^ (9^^ - 2x - 8) Hs4{x)Ho^s4{x) 
lUo Z ( 


108 — ^- ^ ' 2^^ \ / — 

70 70 

— (a; - 6)a;i7o,s4(l)77o,s4(a;) + —(a; - 2)^770 ,s4(1)77o_i,m(1) 


—^{x — 2)^i7o,s4(l)77o ,-1,54(2;)-^ log(2)(a; - 3)(a; + l)i7o,o,s4(l) 

35 35 

- ^ (9a;^ -2x-8) i7o,o,54(l) + ^ (9a;^ - 2x - 8) Hofi,s4{x) 

+ log(2)(a; - 3)(a; + l)i7o,54,54(l) + p (9a;^ - 2x - 8) i7o,*4,54(1) 

- p (9a;^ -2x-8) i7o,*4,54(2;) + ^^(2^ - 3)(2; + l)i7_i,0,0,54(1) 

140 140 

- — (a;-3)(a; + l)i7_ 1,0,0,54(2;) H—^(2^ — 3)(a; + l)i7_i^o,54,54(l) 

y y 

140 140 

-^(2; - 3)(a; + l)i7_i,0,54,54(2;) + -^(2; - 3)(a; + l)i7_i,*4,0,54(1) 

140 140 

- ^(2; — 3)(a; + l)-f7_i^*4_o,s4(2;) H —~ + l)77o ,— 1,54,54 (1) 

140 140 

-^(2; - 3)(a; + 1)770,-1,54,54(2;)-^(2; - 3)(a; + l)77o,o,-i,54(l) 

y y 

140 70 

H—~ 1 ) 770 , 0 ,- 1 , 54 ( 2 ;) —^( 2 ; — 2)^77o,- 1 , 0 , 54 , 54 ( 1 ) 

70 70 

+ “^(^ ~ 2)^77o,- 1 , 0 , 54 , 54 ( 2 ;) —^( 2 ; — 2)^77o,- 1 , 54 , 0 , 54 ( 1 ) 

70 140 

+ ~ 2 )^ 770 ,- 1 , 54 , 0 , 54 ( 2 ;) H —~ 2 )^ 770 , 0 ,- 1 , 0 , 54 ( 1 ) 

140 140 

- —(a; - 2)2770,0,-1,0,54(2;) - -^(^ - 2)2770,0,-1,54,54(1) 

140 140 

H~ 2)2770,0,-1,54,54(2;) H—3“(^ ~ 2)2770,0,0,-1,54(1) 
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i4n 7n 

- —(x - 2 fHo,o,o,-iMi^) + y (^ - 2)2i/oAo,.4,-i(l) 

70 70 

—“ 2)^ifo,o,o,s4,-i(a^) + y ~ 2)^i/o,o,s4,o-i(l) 

70 4 

- y (x - 2)^iJo, 0 ,^ 4 , 0 - 1 ( 2 ^) - y + 1) (4a;^ + 400a; - 1489) C 2 H_i{x) 

140 140 

+ y“(^ ~ 3)(a; + l)i^o,o,^ 4 (l)-f^-i(a;)-^(a; - 3)(a; + l)i^o,s 4 ,s 4 (l)-f^-i(a;) 

1 64 

- ——(x - 1) (52112a;2 - 42806a; - 2393137) HAx) +—(x - 2fHQiix)HQ 1 _i(a;) 

1458 ^ ^ 3 ’ ’ ’ 

8 70 

- - 1) (418a;^ - 2777a; + 32665) C, 2 Hi{x) + -^{x - 2)^iJo,o,M(a;)i?o-i(a;) 

z4o y 

8 592 

+ y (50a;^ - 131a; - 165) C 2 ^^o-i(a^) - y^(^ “ 3)(a; + l)H\{x)Ho_i{x) 

14 70 

- + 1) (280a;2 - 1213a; - 5327) - —{x - 2fHo,oM{^)HoMx) 

140 70 

-y (a; — 3)(a; + l)ifo,s 4 (l)-f^o,-i( 2 ^) + y ~ HQ^si,si{^)HQ-i{x) 

4 

+ + 1) (l04a;^ + 1633a; — 7915) H_i{x)Hq _i{x) 

- y (a; - 3)(a; + l)H^i{x)Ho,i{x) + y (a; - l)(a; + 3)H‘I{x)Ho,i{x) 

+ {21656x^ - 106251a;2 - 38232a; + 710568) Ho,i{x) 

4 256 

- — (2373a;2 - 4638a; - 8990) C 2 i^o,i(a;) + ^{x- 2fHo,i{x)Ho,o,-i{^) 

oi Z( 

- ^(a; + 1) (32a;2 - 5567a; + 18887) H_i{x)Ho,i{x) 

- “ 1) (l4a;^ + 194a; - 355) Hi{x)Ho^i{x) + y (51a;^ - 430) i^o-i(a;)i^o,i(a;) 

4 872 

- y^(^ + 1) (l04a;^ + 1633a; - 7915) iJo-i-i(a;) - “ 2 )\ 2 -f^o- 1 -i(a;) 

1184 592 

+ yy(a: - 3)(a; + l)if_i(a;)iJo,-i,-i(a:) + yy(a; - 2)2ifo,-i(a;)i?o,-i,-i(a:) 

64 2 

+ y (a; - 2)^ifo,i(a;)-f/'o-i-i(a;) + + 1) (32a;^ - 5567a; + 18887) iJo-i,i(a^) 

+ ~ 3)(a; + 1 )H_i{x)Hq_i^i{x) - y (a; - 2 )^i^o-i(a^)^^o-i,i(4 

- ^ (I12a;3 + 2943a;2 - 24948a; + 11629) i/o,o,-i(4 “ 

140 32 

+ -^(a^ - 2 )^i?o,s 4 (l)^^o,o-i(a^) + 7 ^(a^ + l)(4a; - 15)iJ_i(a;)iJo,o-i(a^) 


+ y (a; - l)(25a; + 117)Hi{x)Hofi_i{x) - -{x - 2f Ho_i{x)Hofi-i{x) 
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- — (4112a;^ + 17235a;2 + 8928a; - 11415) Ho,o,i{x) 

16 16 

+ 7^ (83a;^ - 316a; + 68) C 2 ^^o, 0 , 1 ( 2 ^) “ 7 r(^ + ~ 23)i7_i(a;)i7o,o,i(^) 

z ( y 

- ^{x - l)(407a; + 1179)i7i(a;)i7o,o,i(2^) + + ^)H-i{x)Ho^i_i{x) 

4 176 

+ — (529a;^ - 2060a; + 620) Ho^i{x)Hofi^i{x) + -^{x - 2fHo_i{x)Hofl^i{x) 

oi Z( 

2 32 

+ + 1 ) (32a;^ - 5567a; + 18887) i7o,i _i(a;) - y (a; - 2)^i7o- 1 ( 2 ^) 

+ 7 ^ (3600a;3 + 28827a;2 - 50220a; + 145505) Ho,i,i{x) + + ‘2fC2Hoxi{x) 


243 

128 


+ - 3)(a; + l)i7_i(a;)i7o,i,i(2^) + ^(x - ^)ix + 3)i7i(a;)i7o,i,i(^) 

64 2 

- y (a; - 2fHo_i{x)Ho,i,i{x) + — (27a;^ + 1836a; + 812) Ho^i{x)Ho^i^i{x) 

11^4 19^ 

- -^{x - 3)(a; + l)i7o.-i,-i,-i(a:) - y^(a; - 3)(a; + l)i7o,-i,-i,i(^) 

- - 3)(a; + l)i7o- 1 , 1 - 1 ( 2 ^) - - 3)(a; + l)i7o-i,i,i(a;) 

- ^{x + l)(4a; - 15)i7o,o- 1 -i(4 “ ^ (^Sa;^ + 840a; - 718) i7o,o- 1 , 1 ( 2 ^) 

Zi ol 

4 (I 79 a ;2 _ 2740a; + 22) 77o,o,o,-i(a^) + 7 ^ (555a;2 - 5808a; - 11551) i7o,o,o,i(a:) 

^ i 


27 

128 


“ 77 + 340a; - 718) i7o,o,i,-i(a;)- —{x - 3)(x + l)Ilo,i-i-i(x) 

oi o 

- ^(a; - 3)(a; + l)i7o,i-i,i(^) “ ^ (201a;^ + 552a; - 1132) i7o,i,o-i(a^) 

o oi 

+ 4 (2847a;2 - 13308a; - 17566) 77o,i,o,i(a;) - ^(a; - 3)(a; + l)i7o,i,i,-i(a;) 

o i o 

2 ii84 

+ — (I905a;2 - 9008a; - 8556) 7Io,i,i,i(^) “ -^ 0 ,- 1 , 0 ,- 1 ,-i(^) 

32, , 2368, , 

+ “ 3 “^^ “ 2 ) ,- 1 , 1 , 0 ,-i(a^)- 2^4 ~ -^o,o,-i,-i,-i(a^) 

568, / N 128, / N 

+ ~ 2 ) i7o,o,-i,o,-i(a;) H—~ 2 ) 7Io,o,-i,i,i(x) 

H—~ 2)^i7o,o,o,-i,-i(a^) ~ 48(a; — 2)^i7o,o,o,-i,i(a^) 

+ ^(a; - 2)2i7o,o,o,o,-i(a;) - 7 ^ (577a;2 - 1580a; + 1296) i7o,o,o,o,i(a;) 

y 2 < 

i 6 

- 48(a; - 2)^i7o,o,o,i,-i(a^) - ^ (59a;^ - 108a; + 236) i7o,o,i,o,-i(a^) 

- 4 (295a;2 - 2996a; - 316) 77o,o,i,o,i(a;) - ^(a; - 2)2i7o,o,i,i,-i(a;) 


20 



— (x — 2)^ifo,i,-i,o,-i(^) —— 2)^ifo, 1 ,- 1 , 0 , 1 ( 3 ^) 

-— (x^ + 4 ) -f^o, 1 , 0 , 0 ,- 1 ( 3 ^) + ^ (71x^ + 4460x + 1780) 77o,i,o,o,i(3^) 

Z( oi 

64 2 

— (x — 2)^i7o, 1 , 0 , 1 ,- 1 ( 2 ^) ~ ^ (l51x^ + 5788x + 2716) i7o,i,o,i,i(2^) 

64 2 

— (x — 2)^i7o,i,i,o,-i(a^) — ^ (61x^ + 1972x + 948) -f7o,i,i,o,i(a^) 

4832, ^ , 16 , . , , , 

2 ^ + 2) 770 , 1 , 1 , 1 , 1 ( 2 :) + — (six + 354x + 58) ^5 

392 2 

— log(2)(x - 3)(x + 1)C3 - ^ (848 x^ + 819x^ + 102132x - 25745) Cs 

2 35 

— (l47x^ + 2564x + 236) (2(3 — 3 "^^ ~ 2 )^i 7 o,s 4 (l)C 3 

4 32 

-(x + l)(445x - 1329)77 _i(x)C 3 “ " 1)(762: + 249)77i(x)C3 

296 424 

^(x - 2)2i7o,-i(x)C3 + ^(^ + 2)^77 o,i(2:)C3 • 


(32) 


A computer-readable version of this equation can be obtained from 91 . In Eq. (32) 


Til = 5 denotes the number of massless quarks and (n stands for Riemann’s zeta function 
evaluated at n. H^(x) where a only has the elements 0 and ±1 denote HPLs [^. In 
case a contains also s4 the corresponding function refers to the iterated integral with 
square-root element introduced in Eq. (29) of Section]^ In Eq. (32) we observe iterated 
integrals up to weight 5. 


Some of the iterated integrals in Eq. (32), which are evaluated for x = 1, can be trans¬ 
formed to combinations of Riemann zeta functions. However, we prefer to leave 77 q:(1) 


since these terms disappear by construction in case Eq. (32) is evaluated for x = 1 


The square root letter occurring in the result for the topology BT3 has already been 
introduced in Ref. 82 , where it was named The corresponding iterated integrals 


occurred in the context of the calculation of three-loop contributions to massive operator 
matrix elements of Ref. (^. Interestingly, using the substitution x —)■ (1 — x')/x'^, the 


integrals involving in Eq. (|32|) can be brought into the form of cyclotomic polylog¬ 
arithms (cf. 


Ref. 1^) and can thus be represented as Goncharov polylogarithms [9^ 
with the sixth root of unity appearing in the indices, more precisely with the alphabet 
{1,0, (—1)^/^}. Furthermore, all functions without a letter (—1) can be reduced to HPLs 
at the cost of a more complicated argument and an increase of the number of terms. In this 
representation the constants introduced via matching at x = 1 are cyclotomic/multiple 
polylogarithms evaluated at the reciprocal of the golden ratio x' = (a/5 — 1)/2. Neverthe¬ 
less, since the 77...^ 4 ... are by construction real and since their numerical implementation 


is straightforward we decided not to rewrite the expression in Eq. (|32|). 
In Ref, 


46 the second term in the threshold expansion for the N^LO corrections to Higgs 


boson production has been computed. Furthermore, for all contributing partonic channels 
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the exact dependence on x is provided for the coefficients of the leading logarithms in 
log(l — x). In Eq. (32) only (some of) the HPLs are divergent in the limit x —>■ 1 since 
in the iterated integrals involving s4 the letter 1/(1 — x) is absent. After extracting 


the logarithmic divergencies of the HPLs we hnd full agreement with the results given 
in Eqs. (2.26) and (2.27) of for the coefficients of the log^(l — x) and log'^(l — x) 
contribution, respectively]^ 


7 Conclusions 

In this paper we have computed a contribution to the third-order partonic cross section for 
Higgs boson production in gluon fusion, namely the sub-process initiated by two quarks 
with different flavour. The numerical impact of this contribution is small. However, we 
have obtained analytic results retaining the exact dependence on the Higgs boson mass 
and the partonic center-of-mass energy. This constitutes a new result since to date only 
an expansion around the soft limit has been presented in the literature. Our Endings 
constitute an important step towards an exact result of all third-order contributions to 
the Standard Model Higgs boson production. 

In the course of our calculation we have mapped all contributing amplitudes to 17 integral 
families. For each family we have constructed a canonical basis and derived the corre¬ 
sponding system of difierential equations. After evaluating the three- and four-particle 
cut initial conditions the difierential equations could be solved in terms of HPLs in all in¬ 
tegral families except one, which required additional letters in the alphabet of the iterated 
integrals. 
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